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Abstract 

The signals observed at the direct detection experiments DAM A, CoGeNT and 
CRESST could be explained by light WIMPs with sizeable spin-independent cross 
sections with nucleons. The capture and subsequent annihilation of such particles 
in the sun would induce neutrino signals in the GeV range which may be ob- 
served at Super-Kamiokande. We determine the rate of upward stopping muons 
and fully contained events at Super-Kamiokande for various possible WIMP anni- 
hilation channels. This allows us to provide strong constraints on the cross section 
of WIMPs with nucleons. We find that the DAMA and CoGeNT signals are in- 
consistent with standard thermal WIMPs annihilating dominantly into neutrino or 
tau pairs. We also provide limits for spin-dependent WIMP nucleus scattering for 
masses up to 80 GeV. These exclude the DAMA favored region if WIMPs annihilate 
even subdominantly into neutrinos, taus, bottoms or charms. 
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1 Introduction 



Throughout the last decades overwhelming evidence for the existence of dark matter has 
been collected. Despite the precise understanding of the gravitational effects of dark mat- 
ter, its nature still remains a secret. Among the most promising dark matter candidates 
is a particle with weak scale interactions and mass (WIMP). Its abundance from thermal 
production in the early universe would naturally match the dark matter abundance. 

For several years great efforts were made to detect WIMPs directly by their inter- 
action with nuclei. Recently, the three direct detection experiments DAMA [ll[2|, Co- 
GeNT [3] and CRESST |1] have announced signals which can consistently be explained 
by elastic scatterings of rather light WIMPs with masses of a few GeV This interpre- 
tation is now challenged by a low-threshold analysis from the CDMS collaboration which 
excludes the DAMA and a large fraction of the CoGeNT favored region in parameter 
space j6]. Nevertheless, one should be cautious to discard the light dark matter hypothe- 
sis as several concerns about the CDMS analysis have been raised [7] . Furthermore, even 
if the CDMS limit persists light dark matter may still explain at least the CoGeNT and 
CRESST resuhsB 

In any case it seems reasonable to consider independent measures to search for light 
dark matter in collider- and astrophysics. In the past, analysis of the Tevatron and 
LEP data lead to constraints on possible WIMP couplings to quarks, gluons and elec- 
trons [814T0]. Observations of the diffuse gamma ray spectrum [TTHT3] . the antiproton 
spectrum [T3] and the cosmic microwave background |T5] were used to put limits on the 
annihilation cross section of light WIMPs. Most remarkably, the latest data from WMAP 
seem to disfavor a standard thermal WIMP in the GeV range unless it annihilates pre- 
dominantly into muons, taus or neutrinos, or its annihilation cross section is velocity 
suppressed [16]. Neutrino telescopes offer a complementary probe of the light dark mat- 
ter hypothesis. They can detect the neutrino signals from dark matter annihilation and 
are especially sensitive to the neutrino and tau channels. 

In this study we consider WIMP capture and annihilation in the sun including also 
a treatment of velocity suppressed annihilation. We use the observed neutrino spectrum 
at Super-Kamiokande to constrain the WIMP nucleon cross section for various possi- 
ble WIMP annihilation channels. While the Super-Kamiokande collaboration analyzed 
upward through-going muons in the dark matter context |17| . we focus on fully con- 
tained events and upward stopping muons. These event categories are known to be more 
relevant for light WIMPs [TBH2U] but have not been used to constrain WIMP cross sec- 
tions in a model-independent way so far j§ (We also note that an earlier analysis on light 
WIMPs j22] which is also used in j23] contains a serious error and should be taken with 
carejl) 

We obtain strong limits on the spin-independent WIMP nucleon cross section which 

-"^ A precise determination of the parameter space where the CRESST excess can be explained by dark 
matter would require further knowledge especially of the neutron and alpha backgrounds. 

^In |2T] upward stopping muons have been used to constrain specific MSSM neutralinos. 

^ With their formula 14 the authors of [22] calculate the WIMP induced rate of upward going muons in 
Super-Kamiokande. But to derive their limits, they compare it with Super-Kamiokande data on upward 
through-going muons. Thereby, they overlook that upward going muons do not necessarily traverse the 
entire detector volume. Indeed, in the energy range considered most observed muons are created and/or 
stopped inside the detector (see figure 1 in [24]). 
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we relate to direct detection experiments. We find that - for standard astrophysical 
assumptions - a thermal WIMP which annihilates dominantly into tau or neutrino pairs 
is excluded as a possible source of the DAMA and CoGeNT excesses unless the cross 
section is velocity suppressed. But even for pure p-wave annihilation at least the DAMA 
region is ruled out for the two channels. We will comment on models which can account 
for the DAMA and CoGeNT signals and are compatible with the combined data from 
indirect dark matter detection. 

We also provide constraints on the spin-dependent WIMP proton cross section for 
masses up to 80 GeV. Compared to direct detection experiments they are more stringent 
for all considered annihilation channels. 



2 Theoretical Framework 

In this study we constrain the WIMP nucleon cross section by considering WIMP induced 
neutrino signals from the sun. The neutrino spectrum from dark matter annihilation 
depends strongly on the underlying particle physics model. We focus on light WIMPs 
< 80 GeV but otherwise keep our analysis as general as possible. Therefore, we derive 
our limits independently for various possible WIMP annihilation channels. 



2.1 Annihilation Channels 

We consider the annihilation of WIMPs into pairs of Standard Model fermions 

XX bb, cc, Tf, vv . (1) 

For the neutrino channel we assume v = {vT- + v^ + Vf,) /"i-, but - due to neutrino oscillations 
- other combinations would lead to similar limits. Annihilation into photons, lighter 
quarks, gluons, electrons and muons is omitted as the corresponding neutrino signals 
would be neglig ibleH We further take into account 

XX ^ 4r , (2) 

where we assume that the energy is distributed equally between the taus. This chan- 
nel is e.g. accessed if XX annihilates into a pair of new gauge or Higgs bosons which 
subsequently decay into taus. We include it as neutrino telescopes are very sensitive to 
WIMPs annihilating into taus and we want to determine whether constraints can be alle- 
viated if the energy is shared between more taus. The neutrino spectra resulting from ([1]) 
and ([2]) were calculated with PYTHIA Monte Carlo (version 8.1) [25]. We checked that 
they agree well with those given in [26|[27]. only for the tau neutrino production in the 
cc channel there is a discrepancy. This can, however, easily be explained as the spectra 
of |2B | [27 ] were calculated with an older PYTHIA version [2S] which uses an out-dated 
value for the semileptonic branching ratio of Df mesons. 



Annihilation into muons would of course induce a large number of neutrinos. Their energy would, 
however, reside below the detection threshold as muons in the sun - due to their long life-time - are 
stopped before they decay. 
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2.2 Annihilation Cross Section 



Throughout this study we assume that dark matter is produced thermally. In this case the 
annihilation cross section crfj-ei can be related to the dark matter abundance. Expanding 
(JUrei in the relative WIMP velocity av^ei = a + bv'^^^ + 0{vf^i) one finds |29j^ 

^Jg^{TF)Tp{a + ?,TF bj m^) 

where denotes the effective number of relativistic degrees of freedom, Tp the freeze- 
out temperature and the WIMP mass. We set Tp = which is a reasonable 
approximation. In case of s-wave annihilation aVj-^i is determined by the first term in 
the velocity expansion, in case of p-wave annihilation by the second term. The thermal 
WIMP abundance matches the dark matter abundance for 

1.8 • 10"^^ cm^ s~^ 

a = r^=^= (s-wave annihilation) , (4a) 



1.2 ■ 10"^*^ cm^ s"-*^ 

b = r^=^= (p-wave annihilation) , (4b) 



where we take g,, from [30] • We will discuss s- and p-wave annihilation separately in the 
next section. 



3 Dark Matter in the Sun 

WIMPs on their way through the sun may get gravitationally bound if they lose energy 
by scattering with nuclei. The trapped WIMPs can annihilate or escape via evaporation. 
A coherent treatment of the three processes which is valid at low WIMP masses and 
includes velocity suppressed annihilation is rare in the literature. Therefore, we find it 
appropriate to discuss them in some detail. Throughout this study we will use the solar 
model AGSS09 |31] for the composition and properties of the sun. 



3.1 Capture 

A WIMP which scatters off a nucleus at radius r inside the sun gets trapped if its 
velocity after the scattering is lower than the escape velocity ?;esc('")- Let us first discuss 
spin-independent WIMP nucleus scattering. For this the differential contribution of nuclei 
i at radius r to the capture rate takes the form j32j 



E 



i?.max 



dV 2m^^l J u 

^-Ji.min 



^This formula is not applicable to models with resonant annihilation, coannihilations or Sommerfeld 
enhancement. 
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Here PQ^i{r) is the mass density of nuclei i at radius r in the sun and /ij the reduced mass 
of the WIMP nucleus system. The WIMP nucleus cross section o"j is defined as 



^. = ^pA'^, (6) 

where Ai is the mass number of the nucleus i, Cp and yUp the WIMP proton cross sec- 
tion and reduced mass. We assumed equal couplings to protons and neutrons. From 
astrophysics it enters the local dark matter density p-^ ^ 0.3 GeV cm~^ and the velocity 
distribution f{u) of WIMPs outside the gravitational potential of the sun. We use a 
Maxwellian distributioijf] shifted by the velocity of the sun Vq: 



fin) 



with Vq = 220 km s ^. Note that the velocity of an incoming WIMP (before scattering) 
inside the sun is simply given by a/m^ + v'^sd'^)- 



Finally, decoherence effects in the scattering enter via the form factor F{Er), where 
Er is the recoil energy of the nucleus. While F{Eji) = 1 for all relevant energies in the 
case of hydrogen, the suppression can be substantial for heavier nuclei like iron. We use 
Gaussian form factors |33] 

\F{En)\' = e-^^/^' (8) 

with 

i?, = (0.9^^/^ + 0.3) fm, (10) 



where is the mass and Ri the RMS charge radius of the nucleus z. It is known that 
Gaussian form factors are less accurate than e.g. Woods-Saxon form factors especially at 
large Eji. They have, however, the advantage that they can be integrated analytically. 
Further, as only the integrated form factors enter ([5]), the error we make is at the percent 
level [M]. The E^ integration runs from the minimal energy transfer -Eij^min required to 
trap the WIMP to the maximal energy transfer -E/j^max allowed by kinematics. One finds 

1 2n^ 

ER,min = , £R,n,ax = ^ {u^ + vl^{r)) . (11) 

The total capture rate is given by 

Ce = 4vrj:/drr^^, (12) 

* 

where Rq is the solar radius. The sum runs over all types of nuclei in the sun. We consider 
those from hydrogen up to nickel. 



^We neglect that the velocity distribution should be truncated at the local galactic escape velocity 
as this has virtually no impact on capture rates. 
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In the case of spin-dependent scattering the WIMP nucleus cross section does not 
experience the coherent enhancement oc A^. Therefore only capture at hydrogen is rel- 
evant. Its differential contribution to capture rate is given by ([5]) with the form factor 
F{Efi) replaced by 1. 

The capture rates for spin-independent and spin-dependent scattering are shown in 
figure [U 




20 40 60 80 

m j,[GeV] 



Figure 1: Solar capture rate for spin-independent and spin-dependent WIMP nucleus scattering. 



3.2 Annihilation 

WIMPs which are trapped in the sun thermalize and, therefore, follow a Maxwellian 
velocity distribution with temperature T^: 



f{v) = y ^t'^e-(-x^^)/(^^x) . (13) 

It turns out that for all relevant masses the WIMPs reside close to the center of the sun. 
In this region the variations of the solar temperature T^ir) and density pQ{r) are small 
and one can approximate both quantities by constants. A reasonable choice is = Tq (f) 
and pQ{r) = ^©(r) where f is the mean WIMP orbit radius. To derive f we use that the 
WIMP number density scales as [35] 

n^(r) = noe-('"^^('-»/^^ . (14) 

For constant density, the gravitational potential with respect to the solar core reads 
0(r) = 27rp0(f)r^G/3 with G being Newton's constant. By averaging r over the WIMP 
density we arrive at the implicit equation 



7r2GpQ (r)m^ 
which can be solved numerically. 
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If we denote the total WIMP number by A^, their annihilation rate is given by 
Tq = with 

^0 = ]^ / 47rr2nJ(r) {av.^i)^ = | ^ ) (crtVei)© • (16) 

We find that for WIMP masses > 1 GeV the rate is well approximated by 

on o — 0.6 GeV\ , , , , 

= 4.5 ■ 10-^° cm-^ ^^^^^^ j {av,,,)^ . (17) 

If we perform the thermal averaging for an annihilation cross section of the form 0"f rei = 
a + h we find 



(crwrei)© = a + 



6T0(f) 



3.3 Evaporation 

Trapped WIMPs may escape from the sun if they gain energy by scattering with nuclei 
inside the sun This process - called evaporation - is highly sensitive to the 

WIMP mass. The evaporation rate per WIMP roughly scales as 

where Iq ^ 4.7 gyr is the age of the sun and mgvap the evaporation mass. The latter is 
defined as the mass for which Eq = l/t©. Note that if exceeds mevap by a few percent, 
evaporation is totally negligible. Conversely, if falls slightly below mevap virtually all 
trapped WIMPs escape the sun via evaporation. The evaporation rate can be estimated 
as 



8 Cgva.p _ Eggc 



2T0(f 



(20) 



Here E^sc denotes the escape energy at the center of the sun and o"evap the evaporation 
cross section. The latter is given by the total cross section of all nuclei interior to the 
radius ro.95 which is defined by TQ{ro,g^) = 0.95 T0(f). Finally the mean WIMP speed v 
for the thermal distribution ffT3l) reads 



(21) 



Note that the estimate (120|) is only accurate to within a factor of three [37] |3 This 
precision is, however, sufficient for our purposes as it translates into an uncertainty in 



^The error comes from the fact that the WIMP velocity distribution in the sun shows deviations 
from a thermal distribution at high velocities which is not taken into account in (|20p . See discussion 
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the evaporation mass of only 3%. We find that the evaporation mass is well approximated 

by 



mevap = mo + 0.32 GeV log^o [j^zw^) ' (22) 

where niQ = 3.5 GeV (mo = 3.02 GeV) in the case of spin-independent (spin-dependent) 
interactions. 



3.4 Total WIMP Number and Annihilation Signal 

The evolution of the total WIMP number in the sun is described by the differential 
equation 

N = Cq- AqN'^ - EqN (23) 

which takes into account capture, annihilation and evaporation. The general solution to 
this equation is rather unhandy, but in case evaporation can be neglected it takes the 
simple form 



N=^^^^MVC^t) (24) 

and the total annihilation rate at present times reads 

Tq = ^AqN' = tanh^ ( y'C^ t^) . (25) 



In the limit i/C^^A^t© ^ 1 equilibrium between capture and annihilation is reached and 
one finds Fq = Cq/2. For a capture rate of lO^^s"^ this requires (crWrei)© ^ 10~^°cm'^s~^. 
Usually, it is assumed that the leading contribution to the annihilation cross section is 
velocity-independent, i.e. atVei = a with a = O{10~^^ cm^ s~^) (s-wave annihilation). In 
this case equilibrium is safely reached. 

A viable exception is the case where av^ei at freeze-out is dominated by the second 
term in the velocity expansion, i.e. av^ei = a + bv'^^^ with a <^h. Then, one has to take 
into account that the velocity distribution of WIMPs in the sun is different from the 
distribution at freeze-out. For pure p-wave annihilation (a = 0) the thermally averaged 
annihilation cross section in the sun is reduced by a factor TQ{f)/Tp compared to freeze- 
out. In this case capture annihilation equilibrium is typically not reached and F© < Cq/2. 
We decided to derive our limits separately for s- and p-wave annihilation. 

Note, however, that even if a -C 6, the small a may still be sufficient to put capture 
and annihilation into equilibrium in the sun. If this is the case the limits from s-wave 
annihilation apply (with some modifications very close to the evaporation mass). 



4 Neutrino Propagation and Interactions 

Neutrino telescopes like Super-Kamiokande are able to detect muon neutrinos from all 
directions of the sky. In order to determine the incoming muon neutrino signal from dark 
matter annihilation one has to consider neutrino oscillations and interactions. 
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4.1 Neutrino Propagation 



The muon neutrino flux which arrives at the earth from WIMP annihilation in the sun 
may be written as 



(26) 



earth 



where (Iq denotes the distance between earth and sun, (dA/^^^/di?,^^)! earth the propagated 
muon neutrino spectrum per annihilation. We use the density matrix formalism as de- 
scribed in [26l[38l|39] to derive {dN^^/dE^^)\eaTth from the initial neutrino spectra in the 
sun (diV,^,^,,/dE,,,^,J|0. 

The evolution equation used for the propagation from the production point in the 
solar core to the earth reads 



d^ 
dr 



NC 



dg 
dr 



e[H, [H, 0]] 



(27) 



cc 



The diagonal entries of the density matrix g can be identifled with the three neutrino 
flavor eigenstates while off-diagonal entries correspond to superpositions of flavors. The 
terms on the right hand side describe neutrino oscillations, neutral and charged current 
interactions and decoherence effects respectively. For their deflnition see Enii2n]l!l The 
evolution equation has to be solved numerically, the density matrix at r = is given 
by the initial neutrino spectra obtained from PYTHIA. Neutrino interactions will be 
discussed in more detail in the next section. 



4.2 Neutrino Interactions 

Neutrino charged and neutral current interactions with nuclei at high energies {E^ > 
10 GeV) are dominated by deep inelastic scattering. In this regime simple analytic ap- 
proximations for the neutrino cross sections exist (see e.g. Appendix C in [10]) • These 
were continuously used in earlier analyses which mostly dealt with neutrino signals from 
WIMPs of several hundred GeV. 

In this study, however, we consider the annihilation of light WIMPs {m^ < 80 GeV) 
which typically induces neutrinos with energies of a few GeV or even below. For low 
energy interactions, the momentum transfer between neutrinos and nuclei is so small that 
the deep inelastic prescription badly fails. In fact, the neutrino cross section at energies 
-^i^ ^ 3 GeV is dominated by quasi-elastic scattering |1T] and single pion production 
which proceeds through an intermediate nucleon resonance Especially the latter 

makes an analytic treatment of neutrino interactions difficult as a series of resonances 
have to be considered. We therefore decided to simulate neutrino interactions with the 
publically available neutrino event generator NEUGEN 3.5.5 [44] which automatically 
takes into account the discussed processes. A solid treatment of low energy interactions 
is especially important for the calculation of neutrino detection rates, but we also take 
the neutrino cross sections which enter the propagation equation (I27p from NEUGEN. 

^We neglect j^t- regeneration in (dp/dr)|cc sls it is subleading for low energetic neutrinos |21) . 
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5 Neutrino Detection at Super-Kamiokande 



Super-Kamiokande is the most sensitive experiment for neutrino signals in the GeV 
range as they arise from hght dark matter annihilation. Therefore, we will not consider 
other neutrino telescopes. Super-Kamiokande is a water Cherenkov detector of cylindrical 
shape (height: 36.2 m, radius: 16.9 m). It is located in the Kamioka-mine (36° 14' N, 
137° 11' E) and surrounded by rock. The detection rates at Super-Kamiokande depend 
on the energy and the zenith ang 6y of the incoming neutrinos. The double differential 
muon neutrino flux from solar dark matter annihilation can be written as 



dEydcosOy dEy d COS 6*1/ 

The function dP/ dcosOy denotes the differential probability for the neutrino beam to 
arrive with angle cosO^, (averaged over the run time of Super-Kamiokande). It can be 
obtained by tracking the orientation of the Super-Kamiokande detector relative to the 
sun taking into account the orbital and rotational movement of the earth. We employed 
the IDL Astronomy User's Library to calculate dP/dcos6i^, the run time of Super- 
Kamiokande was taken from figure 4.16 in |46) . 

Super-Kamiokande is searching for neutrino induced muons in several event cate- 
gories. We will focus on fully contained muon-like events and upward stopping muons 
(see figure [2]). 



Figure 2: Illustration of event types at Super-Kamiokande: fully contained events (left), upward 
stopping muons (right). 



5.1 Fully Contained Events 

An event is categorized as fully contained muon-like if a particle identified as muon 
is created and stopped within the detector. To eliminate backgrounds, several cuts are 
applied by the Super-Kamiokande collaboration. Events are discarded if the initial vertex 
is less than 2 m away from the detector wall or if the momentum of the muon is below 
a threshold momentum p^^th- The latter is p^,th = 200 MeV for single-ring and p^^th = 
600 MeV for multi-ring events [21] which corresponds to threshold energies of -E^^th = 
226 MeV and -E^^th = 609 MeV respectively. Here single-ring stands for events with 
one observed Cherenkov ring in the detector, multi-ring for those with at least two 

^Meant is the angle between neutrino and cylinder axis of Super-Kamiokande. 
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Cherenkov rings. Multi-ring events are induced if the neutrino interaction gives rise to 
further charged particles (pions) apart from the muon. 

We identify quasi-elastic interactions as single-ring all other types of charged current 
interactions as multi-ring. By this we neglect the finite probability that events involving 
pions are "misidentified" as single-ring. As multi-ring events have a stronger momentum 
cut we thus discard slightly too many events. This, however, only makes our limits more 
conservative. 

The fully contained muon-like event rate can be written as 



mx 1 

R''^ = nn,o [ dE, [ dcos^, . „ "^/^ . [ d^x / dE^ 
' ' dEi, a cos 0^ ' I 



-1 V-2m £^M,th /rfq\ 

1 27r ^ ^ 

f If dcr^^O 

X / dcos6'^^ — / d(Pf,^ — — {x,E^,cos9^{9^^,(pf,^)) . 

J Zn J dt^a Cl cos day 



-1 



Here n^i^o is the number density of water molecules and x the position of the interaction 
vertex. The space integral runs over all points inside the detector volume which are more 
than 2 m away from the detector wall. The function 5^'^ = 1 if the muon stops inside 
the detector, otherwise = 0. Taking the track length of muons in water from table 
11-28 in |17], (5^*^ can be determined by simple geometrical considerations. We further 
denote the zenith angle of the incoming neutrino and outgoing muon by 6i, and 6'^ 
respectively. The angles 6^^,^, (p^y are the zenith and azimuth angles of the muon relative 
to the neutrino, i.e. 



cos 9fj,{9^y, ipf,y) = COS 6y cos 6f,y + sin (p^,y sin 6^ sin O/^y . (30) 

If = if^y = this would imply that neutrino and muon are collinear. 

The differential neutrino charged current cross section at water molecules dcr^^° 
is calculated with NEUGEN (see section 14. 2p . Note that in most analyses performed 
so far, the approximation 6*^ = 6y was used. This would simplify (129|) substantially, 
but - unfortunately - the approximation is not valid for the energy range we consider. 
Therefore, we keep track of the finite angle between neutrino and muon. 

As a cross-check for our approach we used the HKKMS fluxes j38j to calculate 
the (non-oscillated) prediction for fully contained events from atmospheric neutrinos 
in Super-Kamiokande I. We found good agreement with |24) concerning the zenith angle 
distribution of events. We systematically underestimate the number of events by a few 
percent. This is what we expect by our implementation of the event cuts (see above). 



5.2 Upward stopping muons 

Super-Kamiokande can also detect muons which come from neutrino interactions in the 
rock surrounding the detector. To separate neutrino induced muons from cosmic ray 
muons only those traveling in the upward direction are considered. Super-Kamiokande 
distinguishes between muons which traverse the entire detector volume (upward through- 
going muons) and muons which stop inside the detector (upward stopping muons). For 
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both the energy threshold is -E^^th = 1-6 GeV, but the latter are more important for 
constraining light dark matter and we focus on these. Their event rate is given by 

n^x 1 

^stop ^ PRock /"dcos^.— / d\ I dE, 

msi J J dE^dcosOy J J 

} 1 7 da^' ^^^^ 

X / dcos^^^ — / dipf,^ (x, E^, cos9^{e^^, ^^^)) 

J 2n J d-fc^ a cos 

-1 

The mass density of rock pRock can be taken from table IV-6 of jl^. In the same reference 
one can find the energy loss of muons in rock. For simplicity we have assumed that rock 
consists entirely of silicoij^. msi denotes the mass of a silicon nucleus and dcxf^ the 
differential neutrino charged current cross section at silicon. The space integration runs 
over the rock surrounding the detector. Again, the function (5^*°p = 1 if the muon stops 
inside the detector and = otherwise. The Heaviside function in the last line of (13T|) 
takes care that downward going muons are discarded. 

Again we used the HKKMS fluxes to verify that we can reproduce the Monte Carlo 
prediction for upward stopping muons in Super-Kamiokande I with good precision. 



6 Limits from Super-Kamiokande 

Super-Kamiokande performed three runs between 1996 and 2007. The number of ob- 
served muons agrees well with the prediction from atmospheric neutrinos, there is no 
hint for an additional flux of neutrinos from the sun. We will now use the data of Super- 
Kamiokande I, II and III to constrain dark matter annihilation in the sun. 

Super-Kamiokande is able to measure the angle under which a muon is produced. 
At very high energies the muon inherits its direction of flight from the parent neutrino. 
Even for neutrinos in the GeV range the angle between neutrino and muon 6*^,^ is seldom 
larger than 30°. Therefore one can put especially strong constraints on the neutrino flux 
from the sun if one restricts the analysis to muons with < 30° where 6^q is the angle 
between muon and sun-earth direction. (Note that O^q = O^v for neutrinos coming from 
the sun). 

For fully contained muon-like events there exist, unfortunately, only data on the total 
number of events. From figure 1 in [l9] we extract that 8596 events of this type have 
been observed during the time t^'^ = 2806 d. The Monte Carlo prediction for atmospheric 
neutrino events was 86100 This translates into a 90% Poisson upper limit on the (dark 
matter induced) event rate 

Rl'i^ = 13.8 yr-i . (32) 

^°The actual composition of rock is hardly relevant for interaction rates. Standard rock is known to 
consist to equal parts of protons and neutrons |47j . All nuclei Ai with an equal number of protons and 
neutrons have virtually the same ratio ct^' / . . 

^^The Monte Carlo prediction was calculated for the neutrino oscillation parameters set to the best 
fit point of the Super-Kamiokande analysis. A slight variation of 9i3 would, however, only marginally 
affect the prediction as can be seen in the same figure. 
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Figure 3: The Super-Kamiokande 90% upper limits on the spin-independent WIMP nucleon 
cross section for various annihilation channels (equal couplings to proton and neutron are 
assumed). The annihilation cross section is chosen such that the thermal WIMP abundance 
matches the observed dark matter abundance for s-wave annihilation (left panel) and p-wave 
annihilation (right panel). At low WIMP mass the limits arise from fully contained events, at 
higher WIMP mass from upward stopping muons. Below the evaporation mass all constraints 
disappear rapidly. Also shown are confidence regions and limits from direct detection experi- 
ments. 



For upward stopping muons the event distribution with respect to 6*^0 can be taken 
from [SOllST]. Considering only events with 9^q < 30°, the observed number during a 
run time of tS*°P = 2828 d was 53, the Monte Carlo atmospheric neutrino prediction 54 
(see figure 1 of [51]). The corresponding 90% Poisson upper limit on the (dark matter 
induced) event rate is 

^stop,30°^l_24yr-^ (33) 

The prediction for i?^*°P' is given by fl3T|) with the d cos O^y integration running from 
cos(30°) to 1 rather than from —1 to 1. 

The limits (132|) and (133|) translate into limits on the spin-independent and spin- 
dependent WIMP nucleon cross section. The combined limits from fully contained events 
and upward stopping muons are shown together with confidence regions and constraints 
from direct detection experiments in figures [3] and HI We independently treat s-wave and 
p-wave annihilation (see section |3]). 

In the spin-independent case we assumed equal couplings to protons and neutron^, 
in the spin-dependent case couplings to protons onlvP^ For DAMA we calculated the 3(T 
confidence region as described in |52|^ for two different values of the sodium quenching 



^^A different weighting of couplings would merely lead to an overall shift of all limits and confidence 
intervals. This is due to the fact that all nuclei apart from hydrogen have roughly the same ratio of 
protons to neutrons. (And hydrogen does not contribute substantially to the capture in the sun for 
spin-independent scattering.) 

^■^A non-zero coupling to neutrons would not affect the Super-Kamiokande limits, but the limits and 
confidence regions from direct detection experiments. 

^''Spin-dependent scattering is not explicitly covered in but the analysis is analogous. One only 
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Figure 4: Same as figure|3]for spin-dependent scattering. Couplings only to protons are assumed. 



factor gNa. The standard value is gNa = 0.3, but = 0.45 is still within the range of 
experimental uncertainties and would allow to reconcile the DAMA and CoGeNT signals 
(see discussion in j5]). Channeling was assumed to be absent as suggested by |54| . We 
used 600kms~^ for the galactic escape velocity as in the CoGeNT analysis. The CoGeNT 
favored region was taken from the CDMS limit from [B]. We did not try to fit the 
CRESST signal as the collaboration has not yet published their data. 

As can be seen the CDMS limit excludes the DAMA and most part of the Co- 
GeNT favored region]^ However, serious concerns about the CDMS analysis have been 
raised in [7]. If we employ the limits from Super-Kamiokande we find that the CoGeNT 
and DAMA favored regions are excluded for WIMPs which annihilate dominantly into 
neutrino or tau pairs unless the annihilation cross section is velocity suppressed. Even 
for pure p-wave annihilation DAMA remains excluded. The quark channels are less con- 
strained. Our results can be used to restrict possible explanations to CoGeNT and DAMA 
further. Observations from the cosmic microwave background already disfavor ther- 
mal WIMPs with < 8 GeV which annihilate into quarks, gluons, photons or electrons. 
This statement, however, only applies if the annihilation cross sections at freeze-out and 
CMB formation are equal. A possible loop-hole remains velocity suppressed annihilation. 

Summing up the limits from indirect dark matter detection, the annihilation channels 
of light thermal WIMPs are tightly constrained. If we require a simultaneous explanation 
to CoGeNT and DAMA, dominant annihilation into muons remains an option. If the 
annihilation cross section is velocity suppressed all annihilation channels apart from the 
neutrino and tau channels remain viable. In [58] we discussed e.g. a model with p-wave 
annihilation of WIMPs into light Higgs pairs which subsequently decay into taus and 
charms. For this case the 4r-limit alleviated by a factor corresponding to the branching 
fraction into taus (~ 1/2) applies, i.e. the model can explain the CoGeNT and DAMA 
signals without being in conflict with any indirect detection bounds. We also note that 

has to include the spin-structure functions which we take from |53| . 

^'"'Important limits also arise from the XENON 10/100 experiments |551I56| . These depend, however, 
strongly on the low-energy behavior of the scintillation efficiency in liquid XENON which is subject to 
large uncertainties (see discussion in |57|). 
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indirect detection limits may be modified by invoking non-standard astrophysics |59| . 

Let us now turn to spin-dependent WIMP nucleus scattering which was introduced 
to reconcile DAMA with other direct detection experiments [HO]- The DAMA 3a region 
is excluded by PICASSO and COUPP |n2] unless the quenching factor of sodium 
is slightly above the standard value gNa = 0.3. The Super-Kamiokande limits are sub- 
stantially stronger than direct detection limits for all considered annihilation channels. 
Especially the DAMA region is ruled out for WIMPs which annihilate into neutrinos, 
taus, bottoms or charms, even if the branching fraction is only at the few percent level. 

7 Conclusion 

In this study we have analyzed WIMP annihilation in the sun and the corresponding 
neutrino signals. We have determined the expected event rate of upward stopping muons 
and fully contained muon-like events at Super-Kamiokande for various possible WIMP 
annihilation channels. This allowed us to derive strong constraints on the WIMP nucleon 
cross section especially for light WIMPs. We took into account the possibility of velocity 
suppressed annihilation and discussed separately the cases of spin-independent and spin- 
dependent WIMP nucleus scattering. Our results were then related to observations at 
direct dark matter detection experiments. 

For the spin-independent case we were able to exclude WIMPs which annihilate into 
tau or neutrino pairs as a possible source of the signals seen at DAMA and CoGeNT. This 
statement holds unless the annihilation cross section is velocity suppressed, but even then 
the DAMA region remains excluded for these channels. If one takes the existing limits 
from WMAP serious a simultaneous solution to DAMA and CoGeNT with a thermal 
WIMP requires either annihilation dominantly into muonj^ or a velocity suppressed 
annihilation cross section. 

For the spin-dependent case we found that the Super-Kamiokande limits are sub- 
stantially stronger than those from direct detection in the considered mass range (m^ < 
80 GeV). The DAMA favored region is ruled out for annihilation into neutrinos, taus, 
bottoms and charms. The constraints - even for velocity suppressed annihilation - are so 
strong that a solution to DAMA with a (spin-dependent) thermal WIMP seems rather 
difficult. 
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